Mode Bifurcation and Fold Points of Complex Dispersion Curves for the 

Metamaterial Goubau Line 
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In this paper the complex dispersion curves of the four lowest-order transverse magnetic modes 
of a dielectric Goubau line (e > 0, > 0) are compared with those of a dispersive metamaterial 
Goubau line. The vastly different dispersion curve structure for the metamaterial Goubau line is 
characterized by unusual features such as mode bifurcation, complex fold points, both proper and 
improper complex modes, and merging of complex and real modes. 



The Goubau line (G-Line) has been known and studied 
since Sommerfeld and later Goubau considered applica- 
tions using non-radiating surface waves on transmission 
lines L, .2] . While Sommerfeld analyzed a long cylindrical 
metallic wire as the transmission line of interest, Goubau 
realized that adding a dielectric outer sheath to the wire 
reduced the radial extent of the electromagnetic (EM) 
field and thus the dimensions of the associated excitation 
device. Interestingly, while the Sommerfeld wave can ex- 
ist only on a conductor of finite conductivity, the Goubau 
wave can exist even when the inner conductor is assumed 
to have perfect conductivity. The G-Line has been inves- 
tigated since the first part of the twentieth century, and 
its guided modes are well known [l|, 0, 0, 0, . As with all 
open waveguide structures, the G-Line spectrum consists 
of a finite discrete set of guided modes with purely real 
longitudinal propagation constants and an infinite con- 
tinuum of radiation modes. Alsopresent on open loss- 
less structures are leaky waves [a, lZ| characterized by 
discrete complex longitudinal propagation constant solu- 
tions to the dispersion equation, but which are improper 
solutions of Maxwell's equations in that these solutions 
decay longitudinally but do not obey the transverse radi- 
ation condition and thus may only be used in restricted 
regions of space. Improper waves are not considered part 
of an open waveguide spectrum and are often referred to 
as "nonmodal" or "nonspectral" Despite this fact, 
leaky waves have found great usefulness in certain appli- 
cations, particularly those related to leaky wave antennas 
. Some authors have referred to the complex solutions 
of the circular dielectric rod and the standard G-Line 
dispersion equations as leaky modes and have considered 
them on a more or less equal footing with the guided 
modes The leaky waves of even the standard G- 

Line are still not well characterized (only the transverse 
magnetic (TM) solutions have been considered in detail 
but on that structure, all complex leaky wave solu- 
tions of the characteristic equation have EM field com- 
ponents that diverge as the radial coordinate increases to 
infinity and are thus improper. 

In this Letter the G-Line geometry (see Fig. [T|) is 
used with a negative index of refraction dispersive meta- 



FIG. 1: Schematic of the Goubau line geometry. The metal 
core of radius a is a perfect conductor. 



material (NIM) [ll|, [l^ replacing the usual dielectric 
layer. Under these circumstances the metamaterial G- 
Line spectrum consists of guided modes, radiation modes, 
improper complex waves, and proper complex modes. In 
the following we consider only the symmetric transverse 
magnetic (TMo„) solutions of the metamaterial G-Line 
(MM G-line). The transverse electric (TEo„) and hybrid 
modes of the MM G-line will be considered subsequently. 

Using the geometry in Fig. [Tl the characteristic 
equation for the G-Line with perfectly conducting in- 
ner conductor and assuming a 9,z,t dependence of 
^-i(me+fz-uit) ^ can be written in the general form, 
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where 7 = /3 -I- ia is the normalized complex longitudinal 
propagation constant. The longitudinal phase, /?, and 
attenuation a are both normalized by feg = <^/c, the free 
space wavenumber. Considering the m = case, Eq. ([1]) 
reduces to T = for TM modes, where. 



T = €j,Zg = 0, 



(2) 



S - tJL^X/w ~ ix^yg, and g = - H^^\s) / {sH^^\s)) . 
We also define, Z — Jq(w)Yq(v) — Jo(w)Fo(w'); V = 
~Mw)Y^{v) -f .h{v)Yo{w), W = -.h{w)Yo{v) + 
Jq{v)Yi{w), and X ^ Ji{w)Yi{v) - Ji{v)Yi{w), with the 
J's, y's and i/'s the usual Bessel functions of the first 
and second kind, and the Hankel functions of the second 
kind, respectively. The dimensionless wavenumbers are, 
V — k2a, w = k2b, and s = k^b. The transverse propaga- 
tion wave numbers are given by kj = ko\/ ejf^j — 7^, for 
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FIG. 2; Dispersion for standard G-line 

j = 2, 3. For these TMo„ modes, the EM fields in the jth 
region are generaUy written as, Hj = Hj^gOe~^'-''^~'^*\ 
and Ej = {Ej,pp + Ej^zz)e^^^''^~'^^\ for propagation 
in the positive z direction. Maxwell's equations and 
the cylindrical symmetry result in the following trans- 
verse fields, HjQ = —iejko/ K^dEj^z/dp, and Ejp = 
—i'y/KjdEj^z/dp. The forms of ([T]) - ([2]) above are most 
suitable for emphasizing material parameter changes in 
both e and fi of Regions 2 and 3, although they have been 
derived previously in alternate forms [5] . 

We solve Eq. ^ numerically for its complex zeros as a 
function of frequency by employing a recently developed 
complex root finding algorithm. Fig. [2] shows the four 
lowest-order TMo„ (n = 0, 1, 2, 3) modes as functions of 
frequency (in GHz) for the standard G-Line. Part (a) 
shows the normalized longitudinal phase vs. frequency 
while Part (b) shows normalized longitudinal attenuation 
vs. frequency. In Fig. [2] the parameters €2 — 5, p2 — ^, 
a = 5 mm, 6 = 10 mm were used to compare our results 
with those in Ref. In particular, note the lowest-order 
TMqo mode (black curve) that has no cutoff frequency 
and (2 = 0. Except for this mode, our results are con- 
sistent with Ref. y. We have assumed that the outside 
region is air and thus €3 ^ ps — 1. At a given frequency, 
when the imaginary part of a root is zero, that root rep- 
resents a guided or bound solution. When the imaginary 
part is nonzero, that root represents a complex leaky 
wave solution of ([2]). Upon calculating the EM field com- 
ponents associated with certain leaky wave solutions in 
Fig. [21 all discrete complex solutions of the standard G- 
line (with Region 2 index of refraction greater than one) 
are improper [61, |9]. Because the outside region is air, 
cutoff occurs when 7 = 1 or s = 0. When 1 < f3 < y/b, 
and a = 0, the well-known guided modes result. When 
P < 1, a ^ 0, the waves are leaky. Also note that an 
infinite number of higher order leaky waves can occur at 
lower frequencies These leaky waves have high frequency 
cutoffs but no low frequency cutoffs. Generally these 
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FIG. 3: Dispersion for the MM G-line. The curves in (c) and 
(d) give a magnified view of f3 and a, respectively, near the 
bifurcation point (Point 4). 



waves have |a| ^ 1 and thus are not so far found to 
be particularly useful. In Fig. [51 the real and imaginary 
parts of the longitudinal propagation constant are cor- 
related using color. Thus the lowest-order mode (black) 
has zero attenuation and is purely guided [Fig.[2Ka)], the 
second lowest-order mode (magenta) is leaky at the lower 
frequencies with its attenuation magnitude increasing as 
the frequency decreases (see Fig.[2ja) and (b)). For this 
mode as the frequency increases up to cutoff at about 14.9 
GHz, the associated attenuation decreases in magnitude 
until this mode becomes a guided mode once cutoff is 
passed (magenta — > blue). The turquoise blue and red 
curves show the next higher-order modes which for much 
of the frequency range are leaky with \a\ 3> 1. 

If the dielectric layer of Region 2 is replaced by a neg- 
ative index of refraction MM, the TM mode structure is 
vastly different. The MM is characterized by the follow- 
ing dispersion formulas, /i2(,/) = fJ-i+Fp/ifQ-p-iTf), 
and, e2(/) = eH + f^/{f^ - P - iTf). We take the filling 
factor, F, to be 0.25, and the damping factor, F, to be 
zero. The low and high frequency limits are /i; = — 1 and 
th — —5, respectively. The resonance frequencies are, 
/o = 0.75 GHz, = 0.5 GHz, and fp = 0.25 GHz. 

Figure [3l contains the four lowest-order TM modes as- 
suming a metamaterial G-Line with both 62 and /i2 nega- 
tive over the frequency range of interest. Fig. [3l^a) shows 
the normalized longitudinal phase vs. frequency, part (b) 
shows the normalized longitudinal attenuation vs. fre- 
quency, where |a| is small, and also shows the normal- 
ized longitudinal attenuation on a larger scale. In Fig. [31 
£2 w —5 and /12 ~ —1- The radii, a and 6, are the same 
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as for the standard G-line in Fig. [2l Referring to the var- 
ious mode "paths" or "branches" in Fig. [3l note that the 
green curve now has a leaky part (below about 3 GHz) 
and the portion of the green curve with f3 > 1 now has 
a ^ 0. As the frequency decreases, the attenuation peaks 
in magnitude between 4 and 5 GHz (green curve) and 
then starts to decrease. Thus the green path is a wave 
with a complex 7. This path is also improper (until it 
merges with the black curves at much higher frequency) . 
Thus the fundamental guided mode of the standard G- 
line is replaced by a complex mode that is improper. The 
mode represented by the magenta path is also complex 
with a nonzero until point 4 is reached. However this 
path is proper. At point 4 (/ « 10.025 GHz), there is a 
mode bifurcation [see Fig.jS^c) and (d)], as has been ob- 
served in catastrophe and bifurcation theories. Catastro- 
phe theory attempts to study how the qualitative nature 
of solutions of equations depends on the parameters ap- 
pearing in those equations [3, [S] • These techniques are 
invaluable in analyzing the qualitative changes in system 
dynamical characteristics by small perturbations. From 
this standpoint, the dispersion relationship, T = 0, is 
cast as a nonlinear equation in (3, a, and / with param- 
eters ei, £2, Ml 7 ^2, a, and b. Using catastrophe the- 
ory, areas of unusual qualitative interest such as the fold 
point at Point 4 can be determined. Using numerical 
techniques alone to solve the dispersion relation over a 
frequency range is difficult due to its "dense" set of ze- 
ros at each frequency. Point 4 is a complex fold point (a 
point where, because no material losses in the line have 
been assumed as yet, 7, —7, 7* and —7* all meet at 
the same frequency) [11, [13, [3, [l^ . As the frequency 
increases past the fold point, a ^ along both upper 
and lower black paths. The black paths near points 2 
and 3 are therefore both guided modes immediately after 
bifurcation. 

Thus, in Fig. [3] there are an improper complex set of 
roots, a proper complex set of roots, and a mode bi- 
furcation at a fold point in (a) near / « 10 GHz past 
which there are two guided modes (i.e., cf = on both 
black paths after 10 GHz). By considering the two guided 
paths (2 and 3) where cF = 0, we find that the black up- 
per path always has a slope opposite in sign to the lower 
path until the first cutoff frequency, / « 14.9 GHz, and 
thus is backward traveling in this frequency region. At 
cutoff the lower path near point 3 has zero slope and af- 
ter cutoff, the slope becomes positive. Also after cutoff, 
this path becomes complex with |a| <C 1, as shown by 
the blue curve in Fig. [3I^b). The red and turquoise blue 
curves are quite similar to those in Fig. [2] over the given 
frequency range but they, too, exhibit the bifurcation and 
fold point phenomena at higher frequencies None of the 
above situations occur for the standard G-line. 

Certain effects similar to those in Fig. [3] have been seen 
in both grounded metamaterial slab geometries 
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FIG. 4: Normalized power (see text) as a function of fre- 
quency for a few of the dispersion curves in Fig. |31 and using 
the same color scheme. In (a), the radial power flux increases, 
consistent with the attenuation growth for the mode. In (b), 
the longitudinal power flux in the MM region declines rapidly 
with frequency. The modes in (c) and (d) are strictly guided, 
with the EM energy flowing in the MM mainly for those modes 
in the upper branch of the dispersion curve of Fig. |3] 



ometries [H [H [H [2l|, [ll] . In the case of the MM G- 
line, we have a transition from a complex proper mode 
(magneta curve in Fig. [3)^ a)) to a pair of guided modes 
at higher frequencies. We see in both the normalized 
phase and attenuation the characteristic intersection of 
a parabola and a straight hue [Fig. ^c) and (d)]. The 
phase has a straight line to parabola transition character- 
istic of a fold point in catastrophe theory [l^, [13) [3) [l3l ■ 
Note that the attenuation for the MM G-line in Fig.ISfb) 
has changed sign for the red, magenta, and turquoise blue 
curves when compared to the standard G-line. while the 
attenuation shows a parabola to straight line transition. 

The guided and leaky modes of the G-line and MM 
G-line can be described as slow modes and fast modes, 
respectively. The guided modes are slow in that they oc- 
cur in the spectral region y^ei/Is < /3 < y^eTP^- Since 
P > 1, cP/uj > 1 but u}/{c(3) < 1, thus their phase veloc- 
ity is less than the speed of light. In contrast the leaky 
waves are fast waves. Note that with reference to both 
Fig. [2] and Fig. [3] there is a point at the low frequency 
end of each graph where (although we cut the graphs off 
at ^£2/^2) there are complex solutions of the dispersion 
relation that have /3 > ^£2/i2, and a ^ 0. For the stan- 
dard G-line this region is a "forbidden" region. For the 
MM G-line of Fig. [31 /3 > ^£2^*2 is no longer considered 
forbidden, and "superslow modes" can arise [l^ . 
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The total power through the various regions can be 
described by the normahzed power in the radial and 
longitudinal directions, Pju ~ Pj,a/Ptot, for a ^ p,z, 
and j = 2 or 3. The total power, Ptot is the sum 
of power flux, defined as, Ptot = a\Pi.<y\- The 
power flux through the surface of a cylinder of ra- 
dius R and length L encompassing a segment of the 
structure is determined from the spatially integrated 
time-averaged Poynting vector Sj, with components, 
Sj,,{p,z) = {c/87r)e2"^Re[£;j,pi?*,], or S^,p{p,z) = 
-(c/87r)e2"^Re[£;3,z-ff3*e]. Thus, the longitudinal 
power transmitted through the two circular regions is. 



P,. 



-R 



2aL 



lo pdpSj.z{p,0), and P,, 



— 27r pdpSj^z{p,0)- For the surface normal to p, in- 
tegration gives, Pj^p = 27rPe"^(l/Q;) sinh(aL)5'3_p(P, 0). 
We take R/b — 10. In Fig. IH we illustrate the normal- 
ized power as a function of frequency for several of the 
modes shown in Fig. [S] In panel (a), power flows radi- 
ally, with maxima at frequencies correlated to when \a\ 
peaks for that mode (Fig. [31(b)). Although the power 
flow is predominately in the air region, is is evident that 
the counter- directed longitudinal power flux in the MM 
region declines rapidly towards zero with increased fre- 
quency [panel (b)]. This indicates that power emission 
from the waveguide takes place at a slight angle, i3 from 
the z axis over a relatively small bandwidth. As for the 
proper mode spectra, \Pj.z\ = 1/2 for both regions, but 
counter-directed. At the bifurcation point (/ sa 10 GHz), 
the leaky wave transforms into a strictly guided mode at 
higher frequencies, where the two counter propagating 
regions are clearly identified by their sign [Fig. [DJc) and 
(d)]. The lower branch (see Fig. [3]) in the guided mode 
dispersion clearly maps to power flow that is predomi- 
nately in the air region. This is in contrast to modes 
occupying the upper branch in Fig. [31 where again in 
Fig. and (d) we see the power distribution shift from 
both regions equally to reside completely in the MM at 
higher frequencies. This also verifles that the MM G-line 
has the remarkable property that certain guided modes 
can propagate below the first waveguide cutoff and can 
be of the backward as well as of the forward wave type 



In conclusion, dispersion curve behavior for a dielectric 
versus MM G-Line has been shown to be vastly different. 
The MM G-line dispersion is characterized by unusual 
mode bifurcations, complex fold points a la catastrophe 
theory, and both proper and improper complex modes. 
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